Neutrino dispersion in external magnetic fields 
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We calculate the neutrino self-energy operator E(p) in the presence of a magnetic field B. In 
particular, we consider the weak- field limit eB -C m\, where me is the charged- lepton mass cor- 
responding to the neutrino flavor vi, and we consider a "moderate field" m\ <C eB <C m%/. Our 
results differ substantially from the previous literature. For a moderate field, we show that it is cru- 
cial to include the contributions from all Landau levels of the intermediate charged lepton, not just 
the ground state. For the conditions of the early universe where the background medium consists 
of a charge-symmetric plasma, the pure B-field contribution to the neutrino dispersion relation is 
proportional to (eB) 2 and thus comparable to the contribution of the magnetized plasma. 

PACS numbers: 13.15.+g, 14.60.Lm, 95.30.Cq, 97.60.Bw 



I. INTRODUCTION 

The presence of matter or electromagnetic fields mod- 
ifies the dispersion relation of neutrinos in rather subtle 
ways because these elusive particles interact only by the 
weak force. However, Wolfenstein was the first to rec- 
ognize that the feeble matter effect is enough to affect 
neutrino flavor oscillations in dramatic ways because the 
neutrino mass differences are very small [1], with prac- 
tical applications in physics and astrophysics whenever 
neutrino oscillations are important [2]. The presence of 
external fields will lead to additional modifications of the 
neutrino dispersion relation. There is a natural scale for 
the field strength that is required to have a significant 
impact on quantum processes, i.e. the critical value 



B e = ml/em 4.41 x 10 13 G. 



(1) 



Note that we use natural units where h = c = 1 and the 
Lorentz-Heaviside convention where a = e 2 /Ait rs 1/137 
so that e « 0.30 > is the elementary charge, taken to 
be positive. 

There are reasons to expect that fields of such or even 



larger magnitudes can arise in cataclysmic astrophysical 
events such as supernova explosions or coalescing neu- 
tron stars, situations where a gigantic neutrino outflow 
should also be expected. There are two classes of stars, 
i.e. soft gamma-ray repeaters (SGR) [3, 4] and anoma- 
lous x-ray pulsars (AXP) [5, 6] that are believed to be 
remnants of such cataclysms and to be magnetars [7], 
neutron stars with magnetic fields 10 14 -10 15 G. The pos- 
sible existence of even larger fields of order 10 16 -10 17 G is 
subject to debate [8-12]. The early universe between the 
QCD phase transition ( ~ 10~ 5 s) and the nucleosynthe- 
sis epoch ( ~ 10 _2 -10 +2 s) is believed to be yet another 
natural environment where strong magnetic fields and 
large neutrino densities could exist simultaneously [13]. 

The modification of the neutrino dispersion relation 
in a magnetized astrophysical plasma was studied in 
the previous literature [14-17]. In particular, a charge- 
symmetric plasma with m e <C T <C mw and B < T 2 
was considered for the early-universe epoch between the 
QCD phase transition and big-bang nucleosynthesis. Ig- 
noring the neutrino mass, the dispersion relation for the 
electron flavor was found to be [16, 17] 
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where p is the neutrino momentum and <p is the angle 
between B and p. The first term in Eq. (2) is the domi- 
nating pure plasma contribution [18], whereas the second 
term is caused by the common influence of the plasma 
and magnetic field [16]. The third term is of second or- 
der in (eB/T 2 ) <C 1 but was included because of the large 
logarithmic factor ln(T/m e ) ^> 1 [17]. The dispersion re- 
lation of Eq. (2) applies to both v e and D e without sign 



change in any of the terms [19]. 



The B-field induced pure vacuum modification of the 
neutrino dispersion relation was assumed to be negligible 
in these papers. However, recently this contribution was 
calculated for the same early-universe conditions as de- 
scribed above [20, 21]. The dispersion relation obtained 
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in these papers for both v e and v e can be expressed as 
^ = l + V2G F ^sin 2 0e-^/( 2eS ), (3) 

IpI °^ 

where p± is the momentum component perpendicular 
to the -B-field. It is easy to check that this would be 
the dominant -B-field induced contribution by far and 
thus would lead to important consequences for neutrino 
physics in media [22-24] . The importance of the question 
whether the -B-field contribution to the neutrino disper- 
sion relation is dominant or negligible calls for an inde- 
pendent calculation. 

A literature search reveals that calculations of the neu- 
trino dispersion relation in external -B-fields have a long 
history [25-27]. To compare the different results we in- 
troduce the neutrino self-energy operator that is 
defined in terms of the invariant amplitude for the tran- 
sition v — > v by the relation 

M{v^v) = -V{p)V{p)v{p), (4) 

where p is the neutrino four-momentum. Note that we 
use the signature (+, —,—,—) for the four-metric. Within 
the standard model, the general Lorentz structure of S(p) 
in the presence of a magnetic field can be expressed in 
terms of four numerical coefficients a, 6, c, and d as 

= [ a (P7) + b (P7) || + c (P<Pl) + i d (w)] L . ( 5 ) 

where 7 Q are the Dirac matrices in the standard repre- 
sentation and L = — 75) is the left-handed projection 
operator. The Lorentz indices of four-vectors and tensors 
within parentheses are contracted consecutively. For ex- 
ample, (p<py) = p a tp a pj 13 . Further, <p is the dimension- 
less tensor of the electromagnetic field, normalized to the 
external -B-field, whereas (p is its dual, 



<p a /3 = -^£ a p^(p^ ■ (6) 

Finally, in the frame where only an external magnetic 
field B is present, wc take the spatial 3-axis to be di- 
rected along B. Four- vectors with the indices _L and || be- 
long to the Euclidean {1, 2}-subspace and the Minkowski 
{0, 3}-subspace, correspondingly. For example, p± = 
(0,pi,p2,0) and p|| = (po>0, 0,^3). For any four- vectors 
X and Y we use the notation 

(AT) I, = (X <ptpY) = X Y - X 3 Y 3 , 
{XY)± = (XipyY)=X 1 Y 1 +X 2 Y 2 , 
(XY) = [XY)\\-[XY) ± . (7) 

Perturbatively, the matrix element of Eq. (4) corre- 
sponds to the Feynman diagram shown in Fig. 1 where 
double lines denote exact propagators in the external B 
field. Put another way, the self-energy operator corre- 
sponds to this Feynman graph with the external neutrino 



lines truncated. The motivation for our work is that the 
results obtained by different authors at one-loop level do 
not agree with each other. We anticipate our results in 
Table I where we show b, c and d obtained by previous 
authors and from our calculation detailed below. 
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FIG. 1: Feynman diagram for the field-induced contribution 
to the neutrino self-energy operator. Double lines denote ex- 
act propagators of the charged lepton and the W-boson in an 
external B field. The contribution of the unphysical Higgs 
particles can be neglected in the limit mf <C m%/. 

Turning to the interpretation of the coefficients in 
Eq. (5) we note that a does not have an independent 
meaning because for small neutrino energies, E <C mw, 
when a(p) = const., the first term in Eq. (5) is fully 
absorbed by the neutrino wave-function renormalization. 
The coefficient d corresponds to an induced electric dipole 
moment and as such can be non-zero only in the presence 
of the CP-odd field invariant (FF) — 4E • B. Even in 
this case it is strongly suppressed [26] . 

Therefore, in our case of a pure external B-field only 
the coefficients b and c are relevant for neutrino disper- 
sion. For a massless neutrino the Dirac equation in mo- 
mentum space is 

[(P7) - S(p)] v(p) - . (8) 
The dispersion relation follows from 

det [(jry) - b(pj)\\L - c{pipj)L] = . (9) 
For both v and V this implies the same dispersion relation 

j§ = l+ (fc+y)sin 2 0. (10) 

Actually, in a perturbative sense the quantity c 2 is of 
higher order, taking us to the two-loop level. Therefore, 
to lowest order the dispersion relation of massless neutri- 
nos depends only on b. 

For massive neutrinos the situation is more compli- 
cated because a transverse -B-field induces transitions be- 
tween positive- and negative-helicity states by the usual 
spin precession. Put another way, in the standard model 
neutrinos with mass inevitably have magnetic dipole mo- 
ments. In this situation it is not particularly illumi- 
nating to express the effect of the -B-field in terms of a 
modified dispersion relation because an entirely new phe- 
nomenon appears, the mixing of positive- with negative- 
helicity states. Formally one can still proceed as above 
by S(p) — > m + S(p) in the Dirac equation and obtain the 
new dispersion relation. However, the energy eigenstates 
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TABLE I: Coefficients in Eq. (5) for the neutrino self-energy operator in an external B-field. 



Authors 



Rcf. 



Field strength" b x 



V2 7V 2 

G F 



c x 



V2n 2 



d X 



V2n 2 



McKeon (1981) 



[25] 



+3eB 



+2eB 



Erdas & Feldman (1990) [27] Moderate 



Elizalde et al. (2002) 6 
Elizalde et al. (2004) 6 
Our result (2005) c 
Our result (2005) c 



[20] Moderate 



[21] Moderate 



Weak 



Moderate 



(eBf 
3m^ 

eB 



In- 



+ - 



+ 4 



(eBf 
3m^ 

(eBf 
3mfy 



hi 



~~e~B 



3 

+ 4 



+ 2.542 



eB 
2 



-pl/(2eB) 



+ 



4 

3e£ 
4 

3eB 



°B is "weak" for eB <C m| and "moderate" for m'j <C eB <C m^- 
'Neutrino momentum range < <C eB. 
c Neutrino momentum range < p J <C m,^ . 



are no longer the left- and right-handed helicity states 
but rather a superposition that depends on the magnetic 
field orientation relative to p. 

Therefore, for massive neutrinos the effect of B is bet- 
ter illustrated in terms of the equation of motion of a free 
neutrino state in an external homogeneous £>-ficld. In 
this situation neutrino and anti-neutrino states are not 
connected. The Dirac equation implies for the helicity 



amplitudes v± of a massive neutrino 



(ii) 



where E p = (p 2 + to 2 ) 1 / 2 . The dimensionless B-field 
induced mixing matrix is 
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(12) 



where v — \p\/E is the neutrino velocity. The first term 
would represent at v = 1 the energy shift for massless 
neutrinos of Eq. (10) where we have dropped the higher- 
order c 2 term. Of course, positive-helicity (right-handed) 
massless neutrinos do not suffer an energy shift. 

The second line is identical to the effect caused by a 
neutrino magnetic moment 



mc 
~2B 



(13) 



Therefore, within the standard model the coefficient c is 
implied by the well-known result for the neutrino mag- 



netic moment [28, 29] 



3e Gyvh 
8tt 2 V2 



(14) 



Likewise, c can be extracted from Refs. [30, 31], where 
the neutrino transitions vi «-> Vj in an external elec- 
tromagnetic field were investigated. From Eq. (12) 
and Table I one concludes that for to > 10~ 4 eV x 
MeV) (B/B e ) the neutrino energy shift becomes 
essentially non-diagonal. 

If the neutrino mass is of Majorana type, the four neu- 
trino components in the full Dirac equation are not inde- 
pendent and the mass-implied term vanishes — Majorana 
neutrinos do not have a magnetic moment. On the other 
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hand, once neutrino masses are included one can not 
avoid flavor mixing so that the situation becomes yet 
more complicated by the presence of ordinatry flavor os- 
cillations and flavor off-diagonal magnetic spin preces- 
sions caused by magnetic transition moments which ex- 
ist for both Dirac and Majorana neutrinos. We will not 
entangle our discussion with these complications because 
in the standard model the mass-induced effects all hap- 
pen on the external neutrino legs in the Feynman graph 
of Fig. 1. The self-energy operator itself is indepen- 
dent of neutrino masses, at least to one-loop order. 

We begin in Sees. II and III with the technique to 
calculate the neutrino self-energy operator by using the 
charged-lepton and IF-boson propagators in a magnetic 
field. In Sec. IV we calculate the neutrino self-energy 
contribution from the n-th Landau level of the charged- 
lepton propagator in combination with the exact W prop- 
agator. In contrast to previous assumptions [20, 21] we 
find that it is not enough to use only the lowest Landau 
level. In Sec. V we derive explicit results for the neutrino 
self-energy operator in the limiting cases of a "weak field" 
eB <C mf and a "moderate field" < eB < be- 
fore concluding in Sec. VI. 

II. DEFINITION OF S(p) 

The S matrix element for the transition v — > v corre- 
sponds to the Feynman diagram shown in Fig. 1 where 
double lines denote exact propagators in the presence of 
an external magnetic field. A detailed description of the 
calculational techniques for quantum processes in exter- 
nal electromagnetic fields can be found in Ref. [32] . 

For the charged lepton I we consider a negative electric 
charge Qg = — e < 0. The propagator in the presence of 
a constant and uniform magnetic field is translationally 
and gauge non-invariant [33]. It can be expressed as 

S F (x,y) =e i *( x >y"> S(x - y) . (15) 

Here, S(x — y) is the translationally and gauge invariant 
part of the propagator. The translationally and gauge 
non-invariant phase can be defined in terms of an integral 
along an arbitrary contour as 

*(x,y) = -e f d^(£), (16) 

J X 

where K^) = A»{£) + - y) v . 



The corresponding FF-boson propagator can be repre- 
sented in a similar form 

G F pa {x,y)=e^y)G pa {x-y), (17) 

where $>(x, y) is also given by Eq. (16). 

It is useful to consider the Fourier transforms of the 
translationally invariant parts of the propagators 

S{X) = J^S( q )c-^, (18) 

G P AX) = I ^G pa (q)e^ x . (19) 

The sum of the translationally non-invariant phases of 
the lepton and W propagators Eqs. (15) and (17) in the 
loop vanishes, 

*(x,y) + $(y,x)=0. (20) 

This allows one to extract the amplitude M from the S 
matrix element by the standard method, 

M{v^v) = '-^-m^L J ^L S {q)G fja { q -p) 

x^Liy(p), (21) 

where g is the electroweak SU(2) coupling constant of 
the standard model. Comparing this result with Eq. (4) 
one can express the neutrino self-energy operator as 

^{p) = - 1 -^l a LJ aP {p)^L, (22) 

where 

JM = J ^S( q )G 0a (q-p). (23) 
To calculate J a p{p) we need to unravel the propagators. 



III. CHARGED-LEPTON AND W-BOSON 
PROPAGATOR IN A MAGNETIC FIELD 

For the Fourier transform S(q) of the translationally 
invariant part of the lepton propagator Eq. (18) one ob- 
tains in the Fock proper-time formalism [33] 



S(q) 
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JO 



ds 



cos(/3s) 



exp 



\s\m t - 



tan(/3s) 



[(97)|| +me] 



cos(/3s) - sm(/3s) 



cos(/3s) 



>, (24) 



where [3 = eB and mg is the lepton mass. Similarly, the Fourier transform of the translationally invariant part of the 
VF-boson propagator Eq. (19) can be written in Feynman gauge as [27] 



5 



ds 



cos(/3s) 



exp 



9 tan(/3s) 9 



(w) p <t - (w)/kt cos(2/3s) - (^ pCT sin(2/3s) 



(25) 



Manipulations with the exact expressions Eqs. (24) and (25) are extremely cumbersome. Magnetic fields existing in 
Nature probably are always weak compared with the critical field for the W-boson, m^/e ~ 10 24 G. Therefore, the 
W propagator can be expanded in powers of (3 as a small parameter. We find up to second order 



G p(T (q) 



9 P o 



q 2 - m 2 w 



(3 



2 ifp 



(q 2 -m 2 w ) 2 



2q\ 



ip2 9p ^W^ + W^) +A{wU W^^A + °^ ] - 

(26) 

Likewise, the asymptotic expression for the lepton propagator S(q) is realised when the field strength is the smallest 
dimensional parameter, « « m^. In this "weak field approximation" the charged- lepton propagator can be 
expanded as [34] 



S(q) 



(yy) + m e 
q 2 - m 2 f 
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2i 
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+ (D(f3 3 ) . (27) 



One can see from this expansion that the contribution of the region of small virtual momenta q 2 ~ m 2 -C m^y is 
enhanced in each succeeding term. If the propagator is used for a "moderate field," m 2 <C <C m\y, the expansion 
is not applicable and the exact propagator Eq. (24) must be used. 

When the magnetic field is strong enough, B > Be = mj/e, another possibility is to express the charged-lepton 
propagator as an expansion in terms of Landau levels [35] 



%) = E 2_ m ;_ 2n/3 [(97)|| +m t ] 

n=0 



dn{v) - - (7<pj)d! n (v) 



- (yy)_L 2n 



d n (v) 



(28) 



where v = q\j ' (3 and 

d n (v) = (-l) n e- v [L n (2v) - L n _i(2«)] . (29) 

Here, L n (x) are the Laguerre polynomials with the addi- 
tional definition L_i(x) = 0. 

IV. CONTRIBUTION OF THE LEPTON LOW 
LANDAU LEVELS 

As we have already stressed in the introduction, our fi- 
nal result derived in Sec V below strongly disagrees with 
that of Rcfs. [20, 21]. We think that the disagreement 
arises because these authors use only one lowest Lan- 



dau level in the charged-lepton propagator in the case of 
moderate field strengths which they call "strong fields." 
However, the contributions of the next Landau levels can 
be of the same order as the ground-level contribution 
because in the integration over the virtual lepton four- 
momentum in the loop the region q 2 ~ » (i appears 
to be essential. 

To substantiate this point we calculate the contribu- 
tion to the neutrino self-energy operator from the n-th 
charged-lepton Landau level in conjunction with the ex- 
act W-propagator in the limit p^/m^ <C m\y/(3. Sub- 
stituting the exact W^-propagator Eq. (25) and the n-th 
charged-lepton Landau level from Eq. (28) into Eq. (23) 
wc find 



J °P W - J (27r) 4 



(2tt) 4 q\-m\- 2n(3 
ds 



(97)|| 



dn{v) - - (7^7) d' n (v) 



(97U 2n 



d n (v) 



cos(/3s) 



exp 



(3s 



(9-p) 2 



{w)po ~ {w)pa cos(2/3s) - ip pa sin(2/3s) 

(30) 



Terms with even numbers of 7 matrices were omitted because they are removed by the chiral structure of the operator 
Eq. (22). Next we perform a clockwise rotation in the complex plane s = —it and use the identity 



1 



qf, - m 2 - 2n(3 



dr' exp —t' (mj + 2n(3 — q 2 ^j 



(31) 



G 



These manipulations allow us to rewrite the integral Eq. (30) as 



(P) ~ J (2tt 
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) 4 
drdr' 
cosh(/?7 



(<?7) 



(w)p<T 



2 

(W)pa 



- (<?7)±2n 



d„(t>) 



x exp 



cosh(2/3r) + i ^ pCT sinh(2/3r) 
tanh(/3r) 



r' (m 2 + 2n(3 - qfj - r (m^ - (g - p)f) 



/3 



(9-p) 



(32) 



In the integration over d 4 q = d 2 <7|| d 2 gj_, the integrals over d 2 g|| can be easily calculated because they are of Gaussian 



form. As a result we find 
i 
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exp 



xy 
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qj_ exp 
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tanh(ryx) 



'up pa sinh(2?7x) 

{q-p)l |(p7)h 



x + y 



d n (v) 



\ (l'Pl)dn(v) 



d n (v) 



(33) 



where the dimensionless variables x — m^r and y = m^r' have been introduced as well as the parameters r/ = (3/m\ v , 
£ = ■p^jm^ v ~ pj_/m,w and A = raj/m^. From Eq. (33) follows that the essential region of the x variable is x <~ 1 
due to the exponential e~ x . Given the condition 77 -C 1, the argument of the hyperbolic functions is small, r\x <C 1, 
leading to an obvious simplification. One should also take into account the condition q\ ~ /3 caused by the functions 
d n (v), see Eq. (29), containing the exponential e _tI . For a wide range of the numbers n the exponential in the integral 
over d 2 gj_ is simplified, with the only restriction n <C 1/77 = m 2 ^//?: 



exp 



tanh(?7x) 





(<z-p) 2 



~ exp — x 



/Ti- 



ll' 



x exp — x 



2(<zp)_ 
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w 



(34) 



Here, the first exponential is equal to e~^ x . We consider the value p\ to vary in a very wide range, < p\ <C m^/ '/3. 
The second exponential is equal to unity with a good accuracy, because q\ ~ /3 <C to 2 ^ and (gp)j_ <C rafy. With 
these approximations, the integration over d 2 q^ can be easily performed, 

d n (v) 



J d 2 q^d n (v) = 7r/3(2 - 5 nQ ) , J d 2 q A _d n '(v) = -n P8 n0 , Jd 2 q A _(qj)_ 
The investigated integral acquires the form 



= 0. 



4?(p) 



i/3 



16ir 2 m 2 v 



{Pl)\\ 9 P a<2- 



1 - 2 (7^7) 



x da; dy 
Jx + yf 



exp 



y (2nrj + A) 



Taking into account the smallncss of the parameters r\ and A, one finally obtains for n <C m^//? 



4;\p) = 



i(3 



In 1 
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(35) 



(36) 



(37) 



Substituting Eq. (37) into Eq. (22) we finally find the contribution of the n-th Landau level of the lepton propagator 
to the neutrino self-energy operator 

„2 



E M (v) - G F eB m w 
S iP) - V2 2. 2 ~ 2 h 



Pa 



1 



Pi 



'w 



[(2 - 5 n0 ) (p7)|| - <5„ {pifij)] L . 



(38) 



We conclude from Eq. (38) that, contrary to the treatment of Refs. [20, 21], the lowest Landau level does not dominate. 

For higher Landau levels, n > m 2 ^//?, the calculation is more cumbersome. Therefore, using the lepton propagator 
expanded in terms of the Landau levels, with a further summation, is extremely inconvenient. It is much simpler to 
take the exact lepton propagator in the form of Eq. (24). This approach is used in Sec V below. 

In Ref. [21] a numerical test of the lowest Landau level domination was made. However, we believe that the results 
of this test are misleading. In the right-hand side of Fig. 2 of Ref. [21], corresponding to strong fields eB ~ 10 4 m 2 , the 
value yc a defined in Eq. (89) of Ref. [21] tends to unity. From this behavior it was concluded that the lowest Landau 
level domination worked well in this region. However, the left-hand side of that plot, corresponding to weak fields 
eB ~ 10~ 4 m 2 , would have to coincide with, but strongly contradicts, the weak-field result of Ref. [27]. An accurate 
analysis of Eqs. (89) and (90) of Ref. [21] shows that a precise cancellation of the two infinities arises, whereas the 
rest appears to be of order eB/m^ < 10~ 6 for the relevant field values, but not of order unity as claimed there. 
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V. NEW CALCULATION OF S(p) 

We begin our calculation of S(p) with the simpler case of a "weak field" where B defines the smallest energy scale of 
the problem, eB <C to 2 <C rn^y. As an essential simplification we use the field expansions of the Fourier transformed 
W and lepton propagators. Substituting Eqs. (27) and (26) in Eq. (23), and assuming in addition that p]_ <§C m 2 ^, 
we find an expansion of J a p{p) in powers of the field strength, 



16tt 2 



iB 



'■w 



- 2 SW* (P<Plh5 + <Pa0 {Pi) 



+ i 



eB 



71) 



W 



( 2 to 2 3 \ 11 

9a/?(P7)|| (^3 - + i( ^"/3 (P^7)75 - (w)a/3 (P7) ! + •••• (39) 



The dots include terms having the structure g a p (pj) and contain, in particular, the ultraviolet divergence, to be fully 
absorbed by the neutrino wave-function renormalization. They also include terms with an even number of 7 matrices 
that are removed by the chiral structure of 

Substituting Eq. (39) into Eq. (22) we finally find 



S(P) 



V2 47T 2 



3e(pFj) 



e 2 (pFFj) (A 



m 



w 



-In- 



to; 



+ 1 



(40) 



The corresponding coefficients b and c are shown in Table I. 

For a "moderate field," mf <C eB <C to 2 ^ we use the exact expression Eq. (24) for the charged-lepton propagator 
and the expansion Eq. (26) for the W propagator, with the same assumption p\ <C to 2 ^. After a straightforward but 
cumbersome calculation we find 



J a fi{p) 



16tt 2 



eB 



-7; 9c() (W7)75 + Vaf3 (Pi) 



+ i 



eB 



•■w 



<l„ > (P- ( l-> " 



eB 



where 7e = 0.577 ... is the Eulcr constant, and 



In 3 + - 7e - 27 + iip a0 (2^7)75 - (<ptp) a p (Pi) 



(41) 



sinh 2 z 3 + 



-0.055. 



(42) 



The presence of the term (3 2 In (3 with (3 = eB in Eq. (41) shows that an expansion of the lepton propagator in powers 
of (3 as a small parameter is not possible. 

Finally we find 



S(P) 



G F 



V2 47T 2 



3e(pF 7 ) 



e 2 {pFF 1 ) (A 



TO/ 



VI' 



-In 

3 eB 



+ 3.389 



(43) 



Again, the corresponding coefficients b and c are shown in Table I. 
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VI. SUMMARY 

We have calculated the neutrino self-energy operator 
S(p) in a magnetic field at one-loop order. Our results 
for the invariant coefficients of Eq. (5) that character- 
ize S(p) are shown in Table I for weak and moderate 
field strengths and are compared to those of previous au- 
thors. Our results strongly disagree with those of the 
recent Refs. [20, 21] where a large effect was found. For 



moderate fields we have shown that considering only the 
lowest Landau level contribution in the lepton propagator 
is incorrect because the contributions of the next Landau 
levels are of similar magnitude. 

It is instructive to reproduce explicitly the energy of a 
v e or v e in the presence of a CP symmetric plasma and 
the simultaneous presence of a magnetic field. We write 
the energy shift in the form of Eq. (2). For a weak field 
eB <C TOg we find 



E 



= 1 



V2G F 



7n 2 T 4 
15 



T 2 eB 



COS ( 



(eB) 2 



sin 



T 2 , m 
In — - — In — 



'■w 



w 



(44) 



Interestingly, the logarithmic B-field induced plasma 
term in the third term and the logarithmic pure i?-ficld 
term add to ln(T 2 /m 2 v ) so that no electron-mass depen- 
dence remains. For a moderate field m 2 <C eB <C 
the third term is 

+ 8 i n 2 <f>L _ i n Ok _ 3.542) . (45) 
2tt 2 "% \ vn 2 eB ) 

In this case an electron-mass dependence remains. 

In a plasma, the pure B-field term is comparable to the 
logarithmic contribution of the £>-ficld induced plasma 
term derived in Ref. [17]. However, these logarithmic 
terms do not seem to be numerically important relative 
to the term linear in eB. 
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